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Abstract 

Using 2D quantum spin-1/2 model as a concrete example, we stud- 
ied the relation between gapless fermionic excitations (spinons) and 
quantum orders in some spin liquid states. Using winding number, we 
find the projective symmetry group that characterizes the quantum 
order directly determines the pattern of Fermi points in the Brillouin 
zone. Thus quantum orders provide an origin for gapless fermionic 
excitations. 

PACS number: 73.43.Nq, 74.25.-q, 11.15.Ex 

1 Introduction 

Gapless excitations are very rare in nature and in condensed matter sys- 
tems. Therefore, if we see a gapless excitation, we would like to ask why it 
exists. One origin of gapless excitations is spontaneous symmetry breaking, 
which gives Nambu-Goldstone bosons.^, |2[ However, spontaneous symme- 
try breaking is not the only source of gapless excitations. 

Recently, it was pointed out that the internal order of a generic quan- 
tum state cannot be characterized by broken symmetries and the associated 
order parameter. A new concept - quantum order - was introduced to de- 
scribe those internal orders. || |], ||] It turns out that quantum orders provide 
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another origin for gapless excitations. The gapless excitations from quan- 
tum order are not always ordinary bosons, they can be gapless fermions or 
gapless gauge bosons. 

To have a concrete description of quantum order, a mathematical ob- 
ject - projective symmetry group (PSG) - was introduced. 0, Q PSG is a 
projective extension of symmetry group. The PSG description carries more 
information than the symmetry group description since the same symmetry 
group can have many different extensions. Thus two quantum states with 
the same symmetry can have different PSG's which correspond to differ- 
ent quantum orders. The PSG description of quantum order generalize the 
symmetry group description of classical orders. It was by studying many 
examples of quantum ordered states characterized by the PSG that we found 
that gapless excitations can also be originated from and protected by quan- 
tum orders. ||, |H| 

The relation between the gapless gauge fluctuations and PSG has been 
discussed in Ref. [Q, |5|. The relation is simple and straightforward. The 
gauge group for the gapless gauge fluctuations in a quantum state is found to 
be the so called invariant gauge group (IGG). IGG is a (normal) subgroup of 
the PSG such that PSG / IGG gives rise to the symmetry group of the state. 
This result suggests that the gauge fields that govern all the interactions in 
nature may have their origin in quantum order. ||] 

The relation between the gapless fermions and projective symmetry 
groups are more complicated. However, through a study of most general 
fermion-hopping Hamiltonian that is invariant under a PSG, one can show 
that the PSG can indeed protect gapless fermions. Q We would like to stress 
that those gapless fermions are not chiral fermions and the gaplessness is 
not a result of chiral symmetry. Thus quantum orders also provide an origin 
for gapless Dirac fermions. ||] (Note that the model may contain no fermions 
at high energies!) 

In this paper, we are going to study the relation between PSG and 
gapless fermion excitations using an index theorem. We will use the method 
discussed in Ref. || . Through this study, we reveal a more direct connection 
between PSG and gapless fermions. We can see clearly how PSG protects 
gapless fermions. The direct connection between PSG and gapless fermions 
in the spin-1/2 model suggests that the relation between gapless fermions 
and quantum orders may appear in much broader models. 

It was pointed out that PSG alone cannot completely characterize quan- 
tum orders. |[] It is important to find other universal properties of quantum 
states, so that we can use those properties to characterize quantum orders 
in further details. Using PSG and index theorem, we find the distributions 
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of the crystal momenta of gapless fermions to have certain topological in- 
variance. Thus the topology of those distributions can be used to further 
characterize quantum orders. 

As a concrete example, we will concentrate on a quantum spin-1/2 model 
on a 2D square lattice. Despite the absence of any solid proof, we believe 
that, depending on the spin Hamiltonian, the 2D spin-1/2 model can have 
many different ground states corresponding to different spin liquids. Those 
different spin liquids are described by different quantum orders and have 
different gapless fermion excitations. Those gapless fermions are called the 
spinons, which carry spin-1/2 but no charge. In this paper, we will study how 
the gapless spinons are determined by the PSG that describes the quantum 
order. 



2 Winding number and zero eigenvalues 

First let us introduce winding numbers. Consider a 2n dimensional hermi- 
tian matrix H(k) that depends on two parameters k = (k x , k y ). We assume 
that H(k) anticommutes with 7 s , where 




Thus 



has a form 



M < k) s wliF (2) 

M = ( Jj" x " ) (3) 



where | det(/i)| = 1. Now we see that if the phase of det(/i) winds by 2tt 
around a loop in k space, then det(i/(k)) must vanish somewhere inside the 
loop and H has zero eigenvalues at those points. Thus we can determine 
the existence of zero eigenvalue for H by calculating the winding number 
around a loop. To calculate the winding number, we introduce 

a x = — idk x lndet(/i) 

a y = -id ky lndet(h) (4) 



We find 



-^T^M-^M (5) 
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The winding number around a loop C is denned as 

N = {2-kY 1 j dk-a (6) 

which is quantized as an integer. We note that N does not change as we 
deform the loop, as long as the deformation does not pass any point where 
det(i7(k)) = 0. A non-zero winding number around a loop implies the 
existence of zero eigenvalues of H inside the loop. 

We note that, for a H(k), the winding number along a small loop around 
a zero of det(i?(k)) is ±1. Due to the topological property of the winding 
number, a zero of det(i/(k)) cannot just appear or disappear by itself as 
we deform the i?(k), as long as the condition py , fl] = is not violated. 
Zeros can only be created/annihilated in groups with vanishing total winding 
number. 



3 Gapless points in certain spin liquids 

In the slave-boson approach to 2D spin-1/2 system on a square lattice, [0, 
, [llj] a spin liquid state is described by the following mean-field 
Hamiltonian 

Hmean = ' £ + h.C.) + £ ^ (i)aj,rty(i) (?) 

<ij> i 

where we have used the notations in Ref. Here, I = 1,2,3, t 1 are the 
Pauli matrices, and «ij are complex 2x2 matrices that satisfy 

Tr(iiy) = imaginary, Tr(r'«ij) = real (8) 

The last condition ensure the spin liquid described by the ansatz itjj to have 
spin rotation symmetry, g One can see clearly that H mean is invariant under 
the following SU(2) gauge transformation 

if>(i) - G(i)V(i) 
«y - G(i) u ;j Gt(j) (9) 

where G(i) € SU(2). We also note that the ansatz can be viewed as an 
operator which maps a fermion wave function ip(i) to = 
Such an operator will be called the Hamiltonian whose eigenvalues deter- 
mine the fermion spectrum. The gapless fermions correspond to the zero 
eigenvalue of the Hamiltonian. 
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The PSG of an ansatz ity is simply the symmetry group of the ansatz. In 
this paper we are interested in the relation between the PSG and the gapless 
fermions of the ansatz. Using the notation of H , we consider the set of trans- 
formations U that acts on the lattice, consisting of translation, rotation, and 
so forth. This induces a transformation on uy : uy — ► U(u^) = W[/(i)[/(jV 
When we translate, for example, we may get a wy that looks different from 
the we started with, but the two uy's may in fact be related by a gauge 
transformation. For each U we define a gauge transformation Gu{\)eSU{2) 
which induces a transformation on uy : tiy — > G[/(uy) = G[/(i)uyGjy(j). A 
given ansatz has a symmetry if and only if it is invariant under the corre- 
sponding transformation U followed by a proper gauge transformation Gjj: 

GuU( Uij ) = Ui j (10) 

The group generated by the set of all combined transformations G\j\J 
that leave the ansatz invariant is called the PSG. The IGG is defined as 
a subgroup of the PSG generated by Gj, that is, by the Gu for U = I. In 
other words, the IGG consists of pure gauge transformations that leave the 
form of Uij invariant. IGG is a normal subgroup of PSG which allows us to 
define the coset SG = PSG/ IGG. SG is nothing but the symmetry group 
generated by all the Vs. 



3.1 Gapless fermions in Z2A0013 spin liquids 

Now we are ready to study the direct relation between the PSG and gap- 
less fermions. Let us first discuss a simple symmetric ^-linear spin liquid 
characterized by a PSG labeled by Z2A0013.[[§ The Z2A0013 spin liquid is 
described by a large class of ansatz that are invariant under the Z2A0013 
PSG 0. One simple example of those ansatz is given by 
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u i,i+y 



(XT 1 + 7]T 2 ) 
1 



: +X1T 1 
«i,i+2x = X2T 1 + VlT 2 
"i,i+2y = X2T 1 ~ r/iT 2 



a l Q ^ 0, a 2 / = (11) 

More general ansatz can be found in Ref. The ansatz is invariant under 
the following combined transformations {Go, G X T X , G y T y , Gp x P x , Gp y P y , 
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Gp xy P xy , GtT} that act on the fermion wave functions. They are formed by 
a symmetry transformation followed by a gauge transformation. The sym- 
metry transformations are two translations T x : ip(i x ,i y ) — ► ip(i x ~ anci 
T y \ ip(i x ,i y ) -> ip(i x ,iy ~ 1), three parities i^: rp(i x ,iy) -> ip(-i x ,i y ), P y : 
i>(ix,iy) ^(ix^—iy), and tp(i x ,i y ) — ► ip(i y ,i x ), and a time reversal 
transformation T: T(uy) = — uy.j| The corresponding gauge transforma- 
tions are given by 

0,(1) = T°, G„(i)=T°, 

G ft (i) = r°, Gp y (i) = r°, 
G P . v (i) = ir 1 , G T (i) = ir 3 , 
G (i) = -r°, (12) 

where r° is the 2x2 identity matrix. The above defines the Z2A0013 PSG. 
We note that the combined transformation, such as G X T X is a linear operator 
that act on the fermion wave functions ip. From the operator point of view, 
the invariance of the ansatz implies and requires that all the elements in the 
PSG commute with the Hamiltonian uy, except Gt which anticommutes 
with u\y The IGG, a (normal) subgroup of the PSG, is formed by the pure 
gauge transformations in the PSG. We find that IGG is generated by Go 
and IGG = Z%. Therefore, the above ansatz describes a Z2 spin liquids. 
In the momentum space, the mean-field Hamiltonian becomes 

H(\s) = e(k)r 1 + A(k)r 2 (13) 

The spinon spectrum is given by 



£ ± (k) = ±^ e 2 (k)+A 2 (k) (14) 

where 

e = ag + 2x(cos(k x ) + cos(k y )) 

+2xi(cos(/c :c + k y ) + cos(k x - k y )) 
+2x2 (cos (2 A^) + cos(2k y )) 
A = 2rj(— cos(k x ) + cos(fcy)) 

+277i(cos(2fc a: ) - cos(2k y )). (15) 

We see that r 3 anticommutes with -ff(k). We note that Gt = it^ and the 
Z2A0013 ansatz always change sign under the time reversal transformation 
Gt- Thus, if we choose 7 s = r 3 , 7 5 will anticommute with H(\s) for any 
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Z2A0013 ansatz. In this case, we can define a = (a x ,a y ) via Eq. ^J, which 
allow us to calculate topological invariant winding numbers and determine 
the zeros of -ff(k) for Z2A0013 ansatz. 

Now let us consider the symmetry of H (k) as required by the invariance 
of the ansatz under the Z2A0013 PSG. By considering the action of Gp x Tp x 
and Gp y Tp y on the k-space Hamiltonian H(k), we find that the invariance 
of the ansatz under Gp x Tp x leads to 

U Px H(k x ,ky)U px = H(-k X ,ky), 

U Px = r° (16) 
and the invariance of the ansatz under Gp Tp leads to 



U Py H(k x ,k y )U Py = H(k x ,-k 



Vli 



U Py = r u (17) 

The invariance of the ansatz under Gp xy Tp xy gives us a little less trivial 
result 

U Pxy H(k X ,ky)U pxy = H(ky,k X ), 

U Pxy = ir 1 (18) 

We would like to stress that the above results are valid for any Z2A0013 
ansatz. Together with {t 3 ,H} = 0, we find that the Hamiltonian H(h) for 
Z2A0013 ansatz must take the form Eq. Ql3D with 



a{k x , ky) — e( k Xl ky) — ky) — t(.k y , k x ) 

A(k x , ky) = A(~k X , ky) = A(k X , -ky) = ~ A(ky , k X ) (19) 

In fact, Eqs. (|i~3|) and ( |l9|) define the most general Z2A0013 ansatz. 

Now, let us study the symmetries of a. Under Gp x P x and Gp y P y , we 
have 

M(k X ,ky)=M(-k X ,ky), M(k x ,ky) = M(k x ,-ky). (20) 

Thus 

^xij^xi^y) — ^xi^y)-! ®y(k X )ky) — ®y(. ^rr?^/)? 

a x(k X ,ky) = d X (k X ,—ky), CLy(k X , ky) = —Cly(k X , —ky). (^l) 

Since Gt commute with Gp x P x and Gp y P y , the actions of Gp x P x and Gp y P y 
on a are the same as the actions of P x and P y respectively. 
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Next, we consider the action of Gp xy P xy . From the relation between 
H(k) and M(k), we have 

T l M(k x ,k y y = M(k y ,k x ). (22) 

From the definition 

O'X {k x , ky) 

= -ilTV[r 3 M- 1 (k)a fe M(k)] 

= -ibrr[T 3 M- 1 (lt.)d ka T 1 M(k y , k x y\ 

= i^TiyM- l {ky,k X )d kx M{ky,k x )} 

= —a,y{ky, kx). (23) 

We note that P xy changes (a x (k x ,k y ),a y (k x ,k y )) -> (a y (ky,k x ),ax(ky,kx)). 
From the Z2A0013 PSG, we see that Gp xy P xy and Gt (which is 7 5 ) anticom- 

mute. Therefore Gp xy P xy chang eS (CL X (kxi ky), (ly(kxi ky)) > ((ly(ky, k X ) , flx(^) ) ) ) 

due to the 7 5 in the definition of a. The invariance of the ansatz under 
Gp xy Pxy implies 

Q> X \k X ik y ) = CLy{ky,k x ), &y(k x ,ky) = a, x (k yi k x ) (24) 

We see that the symmetry of a is directly related the algebraic property of 
the corresponding PSG. The symmetry of a is summarized in Fig. ||a 

From the symmetry, we see that the winding numbers vanishes for small 
square loops centered along the lines [(0, 0), (±7r, 0)], [(0, 0), (0, ±7r)], [(it, ±7r), (tt, 0)], 
and [(±7r, vr), (0, tt)]. Therefore, a generic Z2A0013 ansatz does not have any 
gap less fermions along those lines. Now consider a loop C formed by the 
two segment C\ and C2 in Fig. [l|a. C\ and C2 are related through P xy . We 
find that J c dk • a = dk. ■ a. Thus, 2(27r) _1 J Ci dk ■ a is quantized as an 
integer, which is the algebraic sum of the numbers of zeros enclosed by the 
loop formed by C\ and C2. Thus, when 2(2-7r) _1 J c dk • a is an odd integer, 
det(il(k)) must vanish at least at one point on the line ((0,0), (tt, tt)). 

Noticing that a zero of det(iT(k)) typically has ±1 winding number, 
when combined with the symmetry of a, we find a typical distribution of 
the zero of det(i?(k)) may look like Fig. |l]b. The pattern of the winding 
number is determined by noting that Gp x P x map a zero at (k x ,k y ) with 
winding number n to a zero at (—k x ,k y ) with winding number — n. Thus 
action of Gp x P x is the same as the action of ordinary parity transformation 



S 




Figure 1: (a) The symmetry of a in the Brillouin zone. The two arrows 
represent a x x and a y y. (b) A typical distribution of the zeros of det(H(k)). 
A solid dot represents a zero with winding number +1 and an empty dot 
represents a zero with winding number —1. The solid- lines and the dash- 
lines (including the ones on zone boundary) are lines of mirror symmetry, 
whose precise meaning is given in Fig. ||. (c) The curves e = (solid lines) 
and A = (dash lines). The direction of the curves is defined such that the 
right-hand-side of the curve has a positive e (or A). 

P x which change the sign of winding number. This is because Gp and Gp x P x 
commute. Similarly, Gp y P y map a zero at (k x , k y ) with winding number n to 
a zero at (k x , —k y ) with winding number —n. However, Gp xy P xy map a zero 
at (k x , k y ) with winding number n to a zero at (k y , k x ) with winding number 
n. It is strange to see that the P xy parity transformation does change the 
sign of winding number. This is because Gp and Gp xy P xy anticommute. 
We would like to stress that the pattern of the winding numbers and the 
positions of the zeros are directly determined from the algebraic property of 
the PSG. 

For the Z2A0013 state, the pattern of zeros can also be determined from 
Eq. (|l9"l). Zeros appear at the crossing point of the e = and A = curves, 
as shown in Fig. |l|c. 

Physically, each zero correspond to a massless two-component Dirac 
fermion in 1+2 dimensions. We see that in general the Z2A0013 state can 
have An two-component Dirac fermions. In particular, Z2A0013 spin liquids 
can have no gapless fermions if is large enough. 

The topological property of the winding number and the symmetry of a 
imply that the topology of the zeros (TOZ) is a universal property of the 
quantum state. (It was shown that the topology of Fermi surfaces determine 
the quantum orders in free fermion systems. || Here the topology of zeros 
is a special case of topology of Fermi surfaces.) TOZ is not changed by any 
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Figure 2: Type 1 transitions in different symmetry environments. Those 
transitions are named type la, type lb, type lc, and type Id. For each 
case the left and the right panels describe the nodes before and after the 
transition respectively. Each case can be "time reverse" of cause. The 
solid-line and dash-line represent line of mirror symmetry. Zeros on the 
two sides of solid-line have opposite winding numbers and zeros on the two 
sides of dash-line have the same winding number. Note, the square does not 
represent a Brillouin zone. 

small deformations of the ansatz, as long as the PSG is not changed. We see 
that the quantum orders in the spin liquid are characterized by (at least) 
two universal properties: PSG and TOZ.Q PSG alone cannot completely 
characterize the quantum orders. For Z2A0013 spin liquids, the TOZ, by 
definition, is a collection of the following data: (a) a pattern of +1 and —1 
zeros along the line ((0, 0), (tt, 7r)), (b) (N+,N^), the number of (+1,-1) 
zeros inside the triangle ((0, 0), (tt, 7r), (0, 7r)) (not including the zeros on the 
sides and corners), (a) is a universal property because the zeros along the 
line ((0, 0), (tt, tt)) cannot move off the line by itself due to the symmetry of 
H(k). 

Now let us consider the transitions between quantum orders with a fixed 
PSG Z2A0013. Those transitions involve a change in TOZ. However, the 
symmetry of the state is not changed. 

First we would like to have a general discussion of the transitions. By 
examining the possible ways in which TOZ can change that are consistent 
with the symmetry of H(k), we find there are several classes of transitions. 
Some of them are described in Figs. ||, ||, and ||. In the type 1 transition, 
a pair of zeros with opposite winding number is created (or annihilated). A 
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Figure 3: Type 2 transitions in different symmetry environments. Those 
transitions are named type 2a and type 2b. 



type 3 






\ 

• 


/ / 

o 






' / 


• — 

\ 






(a) ' 








4: (a) 


Type 


3 transition. 



type 4 


\ 1 




\ i y" 


y 









number of large dots is twice as large as that of small dots. 



type 2 transition involve three zeros, two of them have the same winding 
number and the third one has the opposite winding number. In the transi- 
tion, three zeros combine into one zero. In type 3 transition, two pairs of 
zeros with opposite winding number are created (or annihilated). In type 
4 transition, four zeros carrying a unit of winding numbers changes a zero 
with negative 2 units of winding numbers into a zero with positive 2 units 
of winding numbers. 

Applying the above results to the Z2A0013 state, we find that the pos- 
sible quantum phase transitions (caused by tuning one parameter) can be 
divided, at least, into the follow 5 classes. Because of the symmetry, we only 
need to describe the motion of zeros near a triangle ((0, 0), (jr, 0), (ir, 7r)). 

(A) A type lb transition on ((0,0), (vr,0)) or ((ir, 0), (ir, ir)) sides of the tri- 
angle. 

(B) A type lc transition on on the line ((0,0), (tt,tt)). 

(C) A type 3 transition at k = (0, 0) or (ir, it). 

(D) A type 2a transition on the line ((0,0), (tt,tt)) . 

(E) A type la transition inside the triangle ((0,0), (7r/2,7r/2), (0, tt/2)). 
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3.2 Gapless fermions in Z2Azzl3 spin liquids 



Next we consider the Z2Azzl3 spin liquid state. One simple example of the 
Z2Azzl3 ansatz is given by:pj 



Wi,i+x = XT 1 ~ W 2 

«i,i+y = + W 2 

^i,i+x+y = -XlT 1 

"i,i-x+y = +XIT 1 

^i,i+2x = u i,i+2y = 



1,2.3 
l 



(25) 



The Z2Azzl3 PSG is generated by 



G x (i) = r°, G y (i) = r°, 

G Pa (i) = i(-) ix+iv r 3 , Gp y {i)=i{-r^yT\ 

G Pxy (i) = ir\ G T {\)=ir\ 

G (i) = -r°, (26) 

In the momentum space, we have 

H(k) = e(k)r 1 + A(k)r 2 (27) 

where 



e = 2x(cos(k x ) + cos(ky)) 

+2xi(cos(A; a; + k y ) - cos(k x - k y )) 
A = 2t](- cos(k x ) + cos(k y )) (28) 

We see that the time reversal transformation Gt gives us again 7 s = r 3 
which anticommutes with -£f(k) for any Z2Azzl3 ansatz. We can use a = 
(a x ,a y ) to calculate the topological invariant winding numbers and deter- 
mine the zeros of H(k) for any Z2Azzl3 ansatz. 

Now let us consider the symmetry of H (k) as required by the the Z2Azzl3 
PSG. We find that the invariance of the ansatz under Gp x Tp x leads to 

Up x H{k x ,ky)U Px = H(-k x + 7T,k y + n) : 

U Px = ir 3 (29) 
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and the invariance of the ansatz under Gp y Tp y leads to 



Up y H(k x ,ky)U Py = H(k x + 7r, —k y + ir), 

= ir 3 (30) 

The invariance of the ansatz under G p„„ Tp , gives us 

U Pxy H(k X ,ky)U pxy = H(ky,k X ), 

U Px = ir 1 (31) 

Thus, the Hamiltonian -ff(k) for Z2Azzl3 ansatz must take the form Eq. 
(|3) with 

e(fe x , fey) = -e(-k x + 7T,k y + n) = -e(k x + 7r, -k y + ir) = e(k y , k x ) (32) 
A (A;*, ^) = -A(-fe x + 7r, fe y + 7r) = -A(k x + 7r, -fey + ir) = -A(k y , k x ) 



In fact, Eqs. (f27j) and ( p2[ ) define the most general Z2Azzl3 ansatz. 

Next, let us study the symmetries of a. From the Z2Azzl3 PSG, we see 
that Gp xy P xy and Gt anticommute. The invariance under Gp xy P xy implies 

®x(.k X ,ky) = Qy(ky, k X ), dy(k X , fey) = Q-X^"\j1^x) 

We also see that Gp x P x and Gt commute and Gp x P x changes (k x ,k y ) — ► 
( — fez + vr, fey + 7r). Therefore, the invariance under Gp x P x leads to 

a x{k x ,ky) = — da; ( — fe-j + 7T, fey + 7T ) , 

a y (k x ,k y ) = a y (-k x + iT,k y + it) (34) 
Similarly, the invariance under Gp y P y leads to 

Q>x\kxi ky) = (i x (k x + vr, fey + 7r), 

a y (k x ,k y ) = -a y (k x + ir,-k y + ir) (35) 

The symmetry of a is summarized in Fig. |H|a 

For a loop C formed by the two segment C\ and C2 in Fig. |5|a, we find 
that Jq dk • a = J c dk ■ a and 2(27r) _1 Jq dk ■ a is quantized as an integer. 
Again, when 2(27r)~ 1 J c dk ■ a is an odd integer, det(i?(k)) must vanish at 
least at one point on the line ((0,0), (it, it)). 

A typical distribution of the zeros of det(H(k)) may look like Fig. ||b. 
The pattern of the winding number is determined by noting that two points 
related by Gp x P x or Gp y P y have the opposite winding numbers and two 
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Figure 5: (a) The symmetry of a in the Brillouin zone for the Z2Azzl3 
state. The two arrows represent a x x and a y y. (b) A typical distribution of 
the zeros of det(iT(k)). A solid dot represents a zero with winding number 
+1 and an empty dot represents a zero with winding number —1. The two 
dots marked by "a" are related by the Gp x P x transformation which changes 
(k x ,k y ) — > (— k x + 7r, k y + tt). The two dots marked by "b" are related by 
the Gp y Py transformation which changes (k x , k y ) — > {k x + TT,—k y + tt). (c) 
The curves e = (solid lines) and A = (dash lines). 

points related by Gp xy P xy have the same winding number (see Fig. ^b). 
The quantum orders in the spin liquid are again characterized by (at least) 
two universal properties: PSG and TOZ. For Z2Azzl3 spin liquids, the TOZ 
is again a collection of the following data: (a) a pattern of +1 and —1 zeros 
along the line ((0, 0), (tt, tt)), (b) (N + ,N_), the number of (+1,-1) zeros 
inside the triangle ((0, 0), (tt, tt), (0, tt)) (not including the zeros on the sides 
and corners). 

Now, let us consider the transitions between quantum orders with a 
fixed PSG Z2Azzl3. Again, by examine the possible ways in which TOZ 
can change that are consistent with the symmetry, we find that the quantum 
phase transitions (caused by tuning one parameter) can have, at least, the 
follow 4 classes. Because of the symmetry, we only need to describe the 
motion of zeros near a triangle ((0, 0), (tt, 0), (tt, tt)). 

(A) A type lc transition on the line ((0,0), (-7r,7r)). 

(B) A type 2a transition on the line ((0, 0), (tt, tt)) . 

(C) A type 3 transition at (0, tt) . 

(D) A type la transition inside the triangle ((0,0), (tt/2,tt/2), (0,tt/2)). 
Since the motion of zeros are restricted by the symmetries of H(\t), 

we find that those motion can only change a quantized winding number 
(27r) _1 J c dk. • a by an even integer, where C3 is the loop (0, 0) — > (0, tt) — > 
(tt, tt) -> (0, 0). Thus (2tt)- 1 J c dk ■ a mod 2 is a quantum number which 
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is the same for all the Z2Azzl3 ansatz. A direct calculation shows that 
(2-7t) -1 Jq d\t ■ a mod 2 is equal to 1. Therefore, there must be at least 
one zero on each of the four diagonal lines ((0,0), (±7r, ±7r)). The Z2Azzl3 
spin liquids at least have 4 two-component Dirac fermions. In general it has 
4 + 8n two-component Dirac fermions. See that PSG alone can some times 
guarantee the existence of gap less fermions. 

3.3 Gapless fermions in Z2Bx2(12)n spin liquids 

The third Z2 spin liquid that we are going to study is the Z2Bx2(12)n spin 
liquid. A simple example of this class of ansatz is given by[Q] 

^i,i+x = iXT° + W 2 

a*' 2 ' 3 = 0. (36) 
The Z2Bx2(12)n PSG is generated by 



G x (i) 


= (" 




G y {\) = r°, 


Gpji) 


= <(- 


") T , 


G Py (i) = i(-yyr 2 , 


G Pxy (i) 


= K- 




G T (i) = (-) i -+ i «r° 


Go® 


= (- 


yx+iy T ^ 





where r ab = (r a + r )/v2- The fermion spectrum for ansatz Eq. |36] is 
determined by 

H = -2 X sm(k x )T + 2i 1 cos(k x )T2 

-2 X sm(k y )T 1 + 2r] cos(k y )T 3 (38) 

where k x G (— 7r/2, vr/2), k y G (— 7r/2,7r/2) and 

T = r° ® r 3 ® r°, r 2 = r 2 ® r 3 ® r°, 

r x = r° ® r 1 ® r 3 , r 3 = r 1 ®r 1 ®T 3 . (39) 

H acts on the following eight-component fermion field ^ = V'k+jrx' ^k+iry' ^k+wx+ny) ■ 
We note that H(h) is not periodic in the Brillouin zone, k x G (— ir/2,ir/2) 
and k y G (— vr/2, 7r/2). It satisfies a twisted periodic boundary condition 

B x H(k)B x = #(k + 7rx) 

B y H(k)B y = H(k + iry) 

B 2 = r° ® r 1 ® r° 

B y = r ®^®^ (40) 
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The fermion spectrum can be calculated exactly and its four branches 
take a form ±£'i(k) and zLE^k). The fermion energy vanishes at two iso- 
lated points k = (0, 0), (0, 7r). Near k = the low energy spectrum is given 
by 

E = ±rrV(x 2 + V 2 ) 2 (k 2 x ~ klY + A X ^klkl (41) 

It is interesting to see that the energy does not vanish linearly as k — > 0, 
instead it vanishes like k 2 . 

Now let us consider the symmetry of a generic Z2Bx2(12)n Hamiltonian 
if(k) as required by the Z2Bx2(12)n PSG. The invariance of the ansatz 
under G X T X leads to 

U x H(k)U x = H(k), 

U x = r^r^r 1 (42) 

and the invariance of the ansatz under G y T y leads to 

UyH(k)Uy = H(k), 

U y = t°®t°®t 3 (43) 

Therefore, H(h) has a form 

"<"> " ( o 4 * 4(k> Loo ) (44) 

Under the Gp x P x transformation, we have 

U Px H{k x ,ky)U ] p x = H(-k X ,ky) 

U Px = r 1 ®^®-!- (45) 
Under the Gp y P y transformation, we have 

U Py H(k X ,ky)Up y = H(k X ,~ky) 

U Py = t^t^t 1 (46) 

To calculate the Gp xy P xy transformation, we note that (— ) txlv = (1 + 
_|_ _ (_)»*+h/)/2 and the gauge transformation (-) lxt y changes 

V>k to ^ 

« (V>k + ^k+Trx + V'k+Try - V>k+7r(x+y) ) (47) 



16 



or 

-> 2 ( T ° ® T° <g) T° + T° <g) T 1 <g) T° 

+T° (g — (g (g T 1 ) # k 

(48) 

We also note that P xy switches V'k+Trx an d V'k+Try hi ^k- Thus it changes 

^k x ,k y 5 ® ^ ® ^ ~t~ ^ ® ^ ® t 3 

+r° (g (g + (g T 2 (g r 2 ) * fey)fe;c (49) 

Therefore, the Gp xy P xy transformation changes 

(t° <g> t° <g> t° + t° <g> t 3 <g> t 3 

+ t° (8i r 1 <g r 1 + r° <g r 2 <g> r 2 ) * fey)fe;c 

= ^ (r 12 (g) r 2 (g) r 2 + t 12 (g (g t 3 

+r 12 r 1 (g r° + r 12 0t°® r 1 ) ^ ky ,k x 
= U Pxy V kytkx (50) 

One can check that (Up xy ) 2 = 1. The Gp xy P xy invariance implies that 

U Pxy H(k x ,k y )U Pxy =H(ky,k x ) (51) 

Next, we study the symmetry of a in order to understand gapless fermions 
for a generic Z2Bx2(12)n ansatz. First we consider the transformation GtT. 
We know that u-^ and hence -ff(k) change sign under GtT. Since the ac- 
tion (—)*<*+*» changes — ► t° (g r 1 (g r 1 *k, we find 7^ = r° (g r 1 (g r 1 
anticommutes with H(k) for any Z2Bx2(12)n ansatz. Since U x and U y com- 
mute with the Hamiltonian H(k) we can construct four different 7 5 's that 
anticommute with H(k). The four 7 5 's are given by 

7L = Hif^^hloUTU^ m,n = 0, 1 (52) 

We can use 7^ n to define four a mn and four topologically invariant winding 
numbers. We note that {joo,U x } = 0, {joo,U y } = and {U x ,U y } = 0. 
Thus, 7q , 7i and jf anticommute with at least one of the U x and U y . 
Since U x and U y commute with HQs.) and hence M(k), we can show that 



17 



the winding numbers from 7q , jf and 7^ are always zero. Since U x and 
U y anticommute, we can choose a basis such that U x = r° <8> t° (8> t 1 and 
U y = r° (8> t° tg) r 3 . In this basis -ff(k) and 7^ have a form M4 X 4 (8) r° 
since they commute with both U x and U y . Therefore, the winding numbers 
from 7^ are always even integers. In the following we will only consider the 
71! -winding number. We will set 

7 5 = 7l 5 x =t°0t 2 0t 3 (53) 

and a = a . 

From the Z2Bx2(12)n PSG, we see that GTG x T x G y T y and Gp x P x anti- 
commute. Therefore, 7 s and Up x anticommute and we find that the Gp x P x 
invariance requires 

( ix{k x ,ky S ) = — ( a x ( k x ,ky)), 

a y (k x ,k y ) = -(a y (-k x ,k y )). (54) 

Similarly, we also find that the Gp y P y invariance requires 

&x(kx-,ky) = {<lx(k x , fey)); 

Q"l/(k X ,ky) = —{ — (ly{k X ,—ky)). (^5) 

From Up xy U x U Pxy = U y , U Pxy ^ U Pxy = 7^, and 

Up xy lt n U Pxy = 1 l m {-r n , (56) 
we find that the Gp xy P xy invariance requires 

&x{kxiky) Uy (kyi fer)> 

a y{kxiky) = ci x (ky,k x ). (57) 

Due to the twisted boundary condition Eq. ^TJ, a is not a periodic function 
in the Brillouin zone, k x £ (— 7r/2,7r/2) and k y £ (— 7r/2,7r/2). Since B x and 
By anticommutes with 7 s , we find 

a(k) = -a(k + vrx) , a(k) = -a(k + Try) . (58) 

This allows us to obtain distribution of zeros and winding numbers over an 
extended Brillouin zone. 

From the symmetry of a, we find a typical distribution of the zero of 
det(H(k)) may looks like Fig. ^b. The pattern of the winding number is 
determined by noting that two points related by Gp x P x , Gp y P y , or Gp xy P xy 
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Figure 6: (a) The symmetry of a in the Brillouin zone, k x € (— ir/2,ir/2) 
and k y £ (— tt/2, vr/2), for the Z2Bx2(12)n state. The two arrows represent 
a x x and a y y. (b) A typical distribution of the zeros of det(H(k)). A solid 
dot represents a zero with winding number +2 and an empty dot represents 
a zero with winding number —2. The large dot at k = (0,0) has a winding 
number 4. 

have the same winding number. The quantum orders in the spin liquid 
are again characterized by (at least) two universal properties: PSG and 
TOZ. For Z2Bx2(12)n spin liquids, the TOZ is a collection of the follow- 
ing data: (a) (N + ,N_), the number of (+2,-2) zeros inside the triangle 
((0, 0), (vr/2, 7r/2), (0, 7r/2)) (not including the zeros on the sides and cor- 
ners), (b) a pattern of +2 and —2 zeros on the diagonal line ((0, 0), (7r/2, tt/2)) 

The total sum of the winding numbers in the Brillouin zone is always 
zero for the Z2A0013 and Z2Azzl3 spin liquids, since their -ff(k) is a periodic 
function in the Brillouin zone. For the Z2Bx2(12)n spin liquid, H(k) satisfies 
a twisted boundary condition Eq. ^0[ The total sum of the winding numbers 
Ntoti m this case, can be non-zero. From the symmetry of a, we see that 
Nfot mod 8 is a topological number that is the same for any two ansatz 
that can be continuously deformed into each other via a path that does not 
change the Z2Bx2(12)n PSG. For the Z2Bx2(12)n spin liquid, the N tot mod 
8 is found to be 4. Therefore, det(H(k)) always has a zero at k = (0,0) 
with a winding number 4 + 8n. A generic zero at an arbitrary k has a 
winding number ±2 which gives rise to two two-component Dirac fermions 
with a linear dispersion E(k) oc |k|. The ±4 winding number for the zero 
at k = (0, 0) implies that there are two gapless two-component "Dirac" 
fermions with a quadratic dispersion relation E(k) oc (k) 2 . We see that 
quadratic dispersion relation at k = (0, 0) is protected by the Z2Bx2(12)n 
PSG and is a generic property of the Z2Bx2(12)n spin liquids. 

Now let us consider the transitions between quantum orders with a fixed 
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PSG Z2Bx2(12)n. We find, at least, the following 6 classes of transitions. 

(A) A type lb transition on the side ((0, 0), (vr/2, 0)) or side ((vr/2, 0), (vr/2, vr/2)) 
of the triangle. Note the winding number changes sign if we shift k by vrx/2 

or vry/2. 

(B) A type lc transition on the side ((0,0), (7r/2,7r/2)) of the triangle. 

(C) A type 4 transition at k = (0,0). 

(D) A type 3 transition (in the extended Brillouin zone) at k = (vr/2, vr/2). 
Note the winding number changes sign if we shift k by vrx/2 or vry/2. 

(E) A type 2a transition on the side ((0,0), (7r/2,7r/2)) of the triangle. 

(F) A type la transition inside the triangle ((0,0), (7r/2,0), (vr/2, vr/2)). 



3.4 Gapless fermions in the UlCnOln spin liquids 

Last, we would like to study gapless fermions in the UlCnOln spin liquid (the 
staggered-flux/(i-wave phase of the slave-boson theory [12, 13]). A simple 
ansatz for the UlCnOln state is given by 

4 = o 

«i,i+x = XT 1 + r]T 2 

Ui,i+y = ~ W % (59) 

More general ansatz can be found in Ref. Q. The UlCnOln PSG is gener- 
ated by 

G X (i) = T , Gy(i)=T , 

G Px (i) = r°, G Py (i) = r°, 
G Pxy {\) = ir\ G T {\) = ir\ 
G (i) = e^)^- 3 , (60) 



Since the Hamiltonian uy, G X T X and G y T y commute, the k-space Hamil- 
tonian -f^(k) can be defined which is periodic over the Brillouin zone, k x £ 
(— vr, vr) and k y E (— vr, vr). Since Gt and the Hamiltonian uy anticommute, 
we can choose 7 5 = —iGt and use such a 7 5 to define a. Now we can find 
the symmetry of a through the algebraic property of the UlCnOln PSG. We 
first note that Gt commute with both Gp x P x and Gp y P y . Thus, the Gp x P x 
and Gp y Py invariance of the ansatz implies 

Q'x(.k x ,ky) 0, X (y k X ,ky), Q,y(k X ,ky) Cly ( k X , fcy ) , 

o-x(kx)ky) = a x (k Xl ky), dy{k x ,ky) = a,y{k x , ky). (61) 
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Figure 7: (a) The symmetry of a in the Brillouin zone. The two arrows 
represent a x x and a y y. (b) A typical distribution of the zeros of det(iJ(k)). 
A solid dot represents a zero with winding number +1 and an empty dot 
represents a zero with winding number —1. (c) The curves e = (solid lines) 
and A = (dash lines). 

Since Gp xy P xy anticommutes with Gt, we have, under Gp xy P xy , 

Q>x{kx-iky) = dyiky^kx)-, Ciy{h X ,ky) = a x(ky, fej ) . (^2) 

From [Gq,H] = 0, we find Lq = (— )^+*» r 3 commute with the Hamiltonian 
Uiy Since Lq and Gt commute and the action of Lq increase k by (vr,7r), 
we have 

ky) = Q>x(k x ~\- TT, ky + 7r), (ly(k X) fey) = dy 

The symmetry of a is summarized in Fig. |^a. Comparing with the Z2A0013 
state, here we have an additional symmetry k — > k + 7rx + 7ry caused by Lq. 
The k Hamiltonian H (k) still has the form Eq. [l3|. 

From the symmetry of a, we find that a typical distribution of the zero 
of det(if (k)) may looks like Fig. 0b. The pattern of the winding number is 
determined by noting that two zeros related by Gp x P x or Gp y P y have the 
opposite winding number and two zeros related by Gp xy P xy have the same 
winding number. Two zeros related Lq (a shift of (ir, tt) in k) also have the 
same winding number. 

The quantum orders in the spin liquid are characterized by (at least) two 
universal properties: PSG and TOZ. For UlCnOln spin liquids, the TOZ 
is a collection of the following data: (a) (N + , iV_), the number of (+1, — 1) 
zeros inside the triangle ((0, 0), (ir/2, 7r/2), (0, tt)) (not including the zeros 
on the sides and corners), (b) a pattern of +1 and —1 zeros on the two sides 
of the triangle ((0, 0), (vr/2, vr/2)) and ((0, vr), (vr/2, vr/2)). 
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Let Ntot be the total sum of the winding numbers in the square ((0, 0), (0, tt), (tt, tt), (tt, 0)), 
which is a quarter of the Brillouin zone. Prom the distribution and the re- 
stricted motion of the zeros, we see that Ntot mod 2 is a topological number 
that is the same for any two ansatz that can be continuously deformed into 
each other via a path that does not change the UlCnOln PSG. For the 
UlCnOln spin liquid, the N to t mod 2 is found to be 1. Therefore, det(i?(k)) 
always has a zero at k = (tt/2, tt/2) with a winding number 1 + 2n. We see 
that the zero at (tt/2, tt/2) is protected by the UlCnOln PSG and a generic 
UlCnOln spin liquids have at least four two-component Dirac fermions at 
k = (±7r/2, ±7r/2). 

Now let us consider the transitions between quantum orders with a fixed 
PSG UlCnOln. We find that the possible quantum phase transitions (caused 
by tuning one parameter) can be, at least, divided into the follow 6 classes. 

(A) A type lb transition on the side ((0,0), (0,7r)) of the triangle. 

(B) A type lc transition on the side ((0, 0), (vr/2, vr/2)) or side ((0, tt), (tt/2, vr/2)) 
of the triangle. 

(C) A type 3 transition at k = (0, 0). 

(D) A type 2b transition at k = (7r/2, tt/2) which changes the winding num- 
ber of the zero at (tt/2, tt/2) from 1 to —1 (or from —1 to 1). 

(E) A type 2a transition on the side ((0, 0), (vr/2, vr/2)) or the side ((0, tt), (tt/2, tt/2)) 
of the riangle. 

(F) A type la transition inside the triangle ((0,0), (7r/2, 0), (tt/2, tt/2)). 
The analysis in this subsection also illustrates that the gapless fermions 

and the TOZ are directly determined from the algebra of the PSG. We can 
obtain those results even without writing down the explicit forms of H(k), 
G T , G Px P x , etc.. 

4 Summary 

In this paper, we studied the relation between quantum orders and gapless 
fermion excitations in a quantum spin-1/2 system on a 2D square lattice. By 
studying the symmetry of H (k) as imposed by the PSG and using the wind- 
ing number, we find that PSG directly determine the pattern of Fermi points 
in the Brillouin zone. Thus the gapless fermion excitations (the spinons) are 
originated from and protected by the quantum orders. 
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